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This article reports on a joint theoretical and experimental study of the Pauli quantum-mechanical 
stress tensor Tapix, y) for open two-dimensional chaotic billiards. In the case of a finite current flow 
through the system the interior wave function is expressed as tp — u + iv. With the assumption that 
u and V are Gaussian random fields we derive analytic expressions for the statistical distributions for 
the quantum stress tensor components Tc/j. The Gaussian random field model is tested for a Sinai 
billiard with two opposite leads by analyzing the scattering wave functions obtained numerically from 
the corresponding Schrodinger equation. Two-dimensional quantum billiards may be emulated from 
planar microwave analogues. Hence we report on microwave measurements for an open 2D cavity 
and how the quantum stress tensor analogue is extracted from the recorded electric field. The 
agreement with the theoretical predictions for the distributions for Tai3{x,y) is quite satisfactory 
for small net currents. However, a distinct difference between experiments and theory is observed 
at higher net flow, which could be explained using a Gaussian random field, where the net current 
was taken into account by an additional plane wave with a preferential direction and amplitude. 

PACS numbers: 05.45.Mt, 03.65.-w, 05.60.Gg 
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I. INTRODUCTION 



Chaotic quantum systems have been found to obey re- 
markable universal laws related to, e. g. , energy levels, 
eigenf unctions, transition amplitudes or transport prop- 
erties. These laws are independent of the details of indi- 
vidual systems and depend only on spin and time-reversal 
symmetries. The universality manifests itself in various 
statistical distribution, such as the famous Wigner-Dyson 
distribution for the energy levels in closed systems, the 
Thomas-Porter distribution for wave function intensities, 
wave function form, conductance fluctuations, etc. (for 
overviews, see e.g. [H, H, H, I3l)- Two-dimensional bal- 
listic systems like chaotic quantum billiards (quantum 
dots) have played an important role in the development 
of quantum chaos. These systems are ideal because they 
have clear classical counterparts. Nano-sized planar elec- 
tron billiards may be fabricated from high-mobility semi- 
conductor hetero-structures like gated modulation-doped 
GaAs/AlGaAs and external leads may be attached for 
the injection and collection of charge carriers i. In 
this way one may proceed continuously from completely 
closed systems to open ones. Here we will focus on open 
chaotic systems in which a current flow is induced by ex- 
ternal means. Simulations for open chaotic 2D systems 
have shown, for example, that there is an abundance of 
chaotic states that obey generalized wave function dis- 
tributions that depend on the degree of openness [1, 0] ■ 
There are universal distributions and correlation func- 
tions for nodal points and vortices [1, 11,^3, UH and the 
closely related universal distributions [g, [I4I and correla- 
tion functions for the probability current density [H, . 

In this article we will focus on the Pauli quantum stress 



tensor (QST) for open planar chaotic billiards and its sta- 
tistical properties. As we will see QST supplements pre- 
vious studies of wave function statistics and flow patterns 
in an important way as it probes higher order derivatives 
(irrespective of the chosen gauge) and thereby fine details 
of a wave function. QST was introduced by Pauli ^^,[16} 
already in 1933 but in contrast to the corresponding clas- 
sical entities for electromagnetic fields and fluids [ij] , for 
example, it has remained somewhat esoteric since then. 
On the other hand, studies of stress are in general an 
important part of material science research and, on a 
more fundamental atomistic level, stress originates from 
quantum mechanics. Efflcient computational methods 
based on electronic structure calculations of solids have 
therefore been developed to analyze both kinetic and con- 
figurational contributions to stress [H, [l^ H^j . The re- 
cent advances in nanomechanics also puts more emphasis 
on the quantum- mechanical nature of stress [2l]|. Fur- 
thermore it features in quantum hydrodynamic simula- 
tions of transport properties of different quantum-sized 
semiconductor devices like Resonant Tunneling Devices 
(RTD) and High Electron Mobility Transistors_(HEMT) 
[23] and in atomic physics and chemistry [2^, [2J] . All in 
all, QST is a fundamental concept in quantum mechanics 
that ties up with local forces and the flow of probability 
density. Hence it is natural to extend the previous studies 
of generic statistical distributions for open chaotic quan- 
tum billiards to also include the case of stress. Our choice 
of planar ballistic quantum billiards is favorable in this 
respect as stress is then only of kinetic origin. Moreover, 
the motion in an open high-mobility billiard may ideally 
be viewed as interaction-free because the nominal two- 
dimensional mean free path may exceed the dimensions 
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of the biUiard itself. In this sense we are dealing to a 
good approximation with single-particle behavior. 

There is an ambiguity in the expression for the stress 
tensor because any divergence-free tensor may be added 
without affecting the forces [1^, [1^ . For clarifying our 
definitions and particular choice, we repeat the basic 
steps, albeit elementary, in Pauli's original derivation of 
his QST [Ulill. If V'(x,t) is a solution to the Schodinger 
equation 



(1) 



for a particle with mass m moving in the external poten- 
tial y, the components of the probability current density 
are 



^ ^ dtp ^djj* 

" 2mi \ dxa dxa 



(2) 



Taking the time derivative of ja and using the the right 
hand side of the Schrodinger equation above to substitute 
dijj/dt, Pauli arrived at the expression 
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where Tap is his form of the quantum-mechanical stress 
tensor 



as well as in open billiards. Distribution functions based 
on measurements were obtained for probability densities, 
currents, and vorticities. In addition, vortex pair cor- 
relation functions have be extracted. For all quantities 
studied BE El complete agreement was obtained with 
predictions based on the assumption that wave functions 
in an chaotic billiard may be represented by a random 
superposition of monochromatic plane waves (s^ . 

The layout of the article is the following. In Sec- 
tion |TT] we outline the meaning of QST by referring to 
Madelung's hydrodynamic formulation of quantum me- 
chanics from 1927 [35]. Section [1111 presents the deriva- 
tion of the distribution functions for the components of 
the QST in 2D assuming that the wave function may be 
described in terms of a random Gaussian field and that 
the net current is zero. Although our focus is on 2D the 
results are extended to 3D as well. Section [IV] deals with 
the distribution of the quantum potential that appears 
naturally in the hydrodynamic formulation of quantum 
mechanics. In Section |V] we present numerical simula- 
tions of transport through an open Sinai billiard with 
two opposite leads and a comparison with the analytical 
Gaussian random field model is made. Microwave mea- 
surements are reported in Section IVII and analyzed in 
terms of the quantum stress tensor. A Berry-type wave 
function with directional properties is introduced in the 
same section to analyze the influence of net currents on 
the statistical distributions for Taf3{x,y). 
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In case of planar billiards V may be put equal to zero and 
it is in that form that we will explore Eq. (|4|). The ki- 
netic Pauli QST is sometimes referred to as the quantum- 
mechanical momentum flux density, see e.g. Ref. [20l |. 
^From now on we will simply refer to it as QST. 

There are obvious measurement problems associated 
with QST for a quantum billiard, among them the limited 
spatial resolution presently available (see e.g. Ref. [27j). 
In the case of 2D quantum billiards there is, however, a 
beautiful way out of this dilemma, a way that we will 
follow here. It turns out that single-particle states ip in 
a hard-wall quantum billiard with constant inner poten- 
tial obey the same stationary Helmholtz equation and 
same boundary condition as states in a flat microwave 
resonator do This means that our quantum billiard 
can be emulated from microwave analogues in which the 
perpendicular electric field takes the role of the wave 
function ip. Since the electric field may be measured this 
kind of emulation gives us a unique opportunity to in- 
spect the interior of a quantum billiard experimentally 
[H, [2§, [sO, [m, m, [11] . Using the one-to-one correspon- 
dence between the Poynting vector and the probability 
current density, probability densities and currents have 
been studied in a microwave billiard with a ferrite insert 



II. THE MEANING OF QST 

One of the earliest physical interpretations of the 
Schrodinger equation is due to Madelung who introduced 
the hydrodynamic formulation of quantum mechanics al- 
ready in 1927 t35|]. This is a helpful step to get a more 
intuitive understanding in classical terms of, for exam- 
ple, quantum-mechanical probability densities and the 
meaning of quantum stress (see e.g. Refs. [H, [13, [ssf). 
Madelung obtained the QM hydrodynamic formulation 
by rewriting the wave function ip in polar form as 



ip(y.,t) = R{x,t) e 



(5) 



The probability density is then p — R^. By introducing 
the velocity v — V5(x, t)/m the probability density cur- 
rent or probability flow is simply j — pv. Intuitively this 
is quite appealing. Inserting the polar form in the Pauli 
expression for Tap in Eq. ([4]) we then have 



^"'^ - A^ 



d'^p 1 dp dp 



dXaXp p dXa dxp 



+ pmvavp. (6) 



There are two qualitatively different terms in Eq. ([B]), a 
quantum-mechanical term Tap that contains the factor h 
and therefore vanishes in the classical limit h ^ 0, plus 
the "classical" contribution pmVaVp which remains in the 
classical limit. Using the notations above Eq. ^ gives 
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the quantum hydrodynamic analogue of the familiar clas- 
sical Navier-Stokes equation for the flow of momentum 
density mp\ 



dpva 



(7) 



Alternatively the Schrodinger equation may be rewritten 
as the two familiar hydrodynamic equations in the Euler 
frame [H, [13, [11 



— + [v • V]v = f/m + F/m, 



(8) 
(9) 



where the external force is due to external potential 

f = -VV, (10) 
and the internal force is due to the quantum potential 



F — — VVqm, Vqm — — 



2m R 



(11) 



Then the internal force can be expressed by a stress ten- 
sor for the probability fluid as 



(12) 



Thus we are dealing with a "probability fluid" in which 
flowlines and vorticity patterns are closely related to 
QST. 



(A' + 1) u{x' ^ y') = and similarly for v{x\ y'). The size 
of the cavity scales accordingly a,s A A' . 

If the shape of the cavity is chaotic we may assume 
that u and v are to a good app roximation random Gaus- 
sian functions (RGFs) [1, Ull with (u^ + w^) = 1 + e^, 
(v"^) — e^(u^), {uv) = and {u) = [v) ~ Q. If u and v 
were correlated we can apply a phase transformation Q 
which makes these functions uncorrelated. Here we use 
the definition 
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(14) 



In what follows we thus use dimensionless derivatives in 
x' and express the QST components in units of the energy 
h^k^/2m. From Eq. dropping the prime (') in the 
expressions from now on, we then have 



and 



^u— — V 



du du dv dv 



dxdy dxdy dx dy dx dy 



(16) 



Two-dimensional case. Let us first consider the distri- 
butions of the stress tensor for a two-dimensional com- 
plex RGF ^j). In the following derivation we assume that 
the net current from one lead to the other is so small 
that in practice we are dealing with isotropic RGFs. We 
therefore have 



(UUx 



1 

'2' 



{UxU:cx) = 0, {ul^) 
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0, 



(17) 



III. DISTRIBUTION OF QST FOR A 
QUANTUM BILLIARD 

We now return to the full expression for the stress 
tensor T^p in Eq. Consider a flat two-dimensional 
ballistic cavity (quantum dot) with hard walls. Within 
the cavity we therefore have = and the correspond- 
ing Schrodinger equation is (A 4- fc^)'0(x, y) = with 
fc^ = 2mE/h^, where k is the wave number at energy 
E. In this case the wave function may be chosen to be 
real if the system is closed and, as a consequence, there 
is no interior probability density flow. The wave function 
normalizes to one over the area A of the cavity. On the 
other hand, if the system is open, for example by attach- 
ing external leads, and there is a net transport, the wave 
function must be chosen complex. Thus 



u + iv, 



(13) 



in which u and v independently obey the stationary 
Schrodinger equation for the open system. In the fol- 
lowing discussion it is convenient to make a substitu- 
tion to dimensionless variables, fcx — > x'. Hence we have 



for the two-dimensional case. The corresponding expres- 
sions for V follow simply by replacing u,Ux,Uxx etc. by 
v/e,Vx/e,Vxx/^ and so on. 

For the component T^x in Eq. (fTSl) we need the follow- 
ing joint distribution of two RGFs [40| 



1 



: exp 



where — {u,v,Ux,Vx,Uxx,Vxx), and the matrix K 



(18) 



). For an isotropic RGF there are only correlations 
{uUxx) , {vvxx)- Therefore the only nontrivial block of the 
total matrix K is the matrix 



1 -1/2 
-1/2 3/8 



K- 



3 4 

4 8 



(19) 



for the RGFs and the matrix Ky — eK^ for the 

two RGFs for v and Vxx- Correspondingly we obtain from 
Eq. (HT 



V8 r 3u2 

/ (m, Uxx) ^ exp \ 

zvr 



8uUj 



8ui 



(20) 
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and 



3^;^ + Svv^a, + 8vl 
2e2 



(21) 



The characteristic function of the stress tensor compo- 
nent T^x is 

e{a) = (e"^-) (22) 

and takes the following explicit form 

0{a) = 8[(1 - m)(l - iea) 

(a - i{V24: + 4))(ea - i(\/24 + 4)) 

(a + i{V24 - 4))(ea - z(\/24 - A))]-^^\{23) 

As a result we obtain for the distribution function 



P{Txx) 



1 

2^ 



6l(a)e" 



-iaT^, 



da. 



(24) 



For e 9^ 1 this integral may be calculated numerically. 
However for e = 1 it might be evaluated analytically. In 
particular, for T^x > we obtain 



P{Txx) = 

and for Txx < 



2 g-(v^-4)T,, 

VG (5-\/24) 



P{Txx) = 



2 g(V24+4)T,, 

7i (5 + V24) 



(25) 



(26) 



The distribution ([26|) is shown in Fig. [T] together with 
results for different e-values obtained by numerical eval- 
uation of the integral p4|) . Note that the distributions 
are here given in terms of (T^x) — 1 + e^- 

To repeat the calculations for the component Txy we 
need the following correlators 

(uUxy) = 0, {UxUxy) = 0, {UyUxy) = 0, {uly) = ^ (27) 

for the 2D case. The correlation matrix turns out to be 
diagonal. Then the characteristic function 

6{a) = 2[(2+(a/2)2)(2 + (ea/2)2)(l + (a/2)2) 

(l + (6a/2)2)]-V2 (28) 

defines the distribution P{Txy). For e = 1 the inte- 
gral ([M)) may, as above, be performed analytically to 
give 



P{Txy) = 2e-2|^-«l - \/2e^2^l^-''l. 



(29) 



The distributions P{Txy) in ((29)l are shown in Fig.[2]for 
the two cases e = and e = 1. Only two cases are shown 
because of the small differences in PiTxy) for different e- 
values. The distributions are in this case given in terms 
of ^(iJJ, where (T^) - |(1 + e^). 
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FIG. 1; (color online) The distribution P{T^^) for e = 1 
(dash-dotted line), e = 0.5 (dashed line), and e = (solid 
line). The stress tensor component T^x is measured in terms 
of the mean value {T^x)- 




FIG. 2: (color online) The distribution P(T^h) for e = 1 
(dash-dotted line) and e = (solid line). The stress tensor 
component T^y is measured in terms of mean value \/ {Tiy). 



Three-dimensional case. In this case the expressions 
in (fT7|) are to be replaced by 



(uUxx) = {uD = ^, 



1 



(uUx) = 0, (UxUxx) = 0, (u^x) = 7 (30) 



and dill) by 

(uUxy) = 0, (UxUxy) = 0, (UyUxy) 



0, {Uly) 



15 



(31) 
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Accordingly the correlation matrix ((T9)) is 
1 



Ku = 



-1/3 \ 9 15 , 

1/3 1/5 i ' " " 4 I 15 45 ' ■ 



The joint probability function of two RGFs u and Ux 
then takes the following form 



45 r + 30mm^^ + 45ul^ 

exp ^ f . [oo) 



2tt 



The characteristic function defining the distribution 

is 



9ia) 



45 



(3/2 - ia)iia + 15/4 + 9^5/4) 
1 

(ia + 15/4- 9\/5/4) 



and, correspondingly, 
P{Txx) ^ 



9>/5-15 r 



15 



(7^5-15) 



for T^t; > 0, and 



(7%/5 + 15) 



9n/5 + 15 rp 

e 4 ^ 



(34) 



(35) 



(36) 



for Ta;a; < 0. Identical expressions hold for the two other 
diagonal components. 

In a similar way we obtain the distribution function 
for the off-diagonal components a ^ (3. For the specific 
case e = 1 we have according to Eq. ([3T|) 



and 



e{a) 



P{Txy) 



[3+(a/2)2][l + (a/2)^ 
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3 1 Tx. y I 



3Vll 



-VT5|Txa| 



(37) 



(38) 



The expression for the other off-diagonal components are, 
of course, identical. 



IV. DISTRIBUTION OF QUANTUM 
POTENTIAL 



The quantum or internal force in Eq. (jlip in the hydro- 
dynamic formulation is defined by the quantum potential 
Vqm- In terms of the RGFs u, v it may be written as 



Vr 



QM 



a' 



(39) 



Vx = 



UUx + VVx 



UUyy + VVyy + + 



-'vv 



UUy + VVy 



The second derivatives might be eliminated using the 
Schrodinger equations for u and v, i.e., Uxx 



^xx ~t- Vyy 

Vqm = 1 
which implies 



- —V. As a result we have 

{uVx ~ VUxY + {uVy - VUyY 



'■vv 



(40) 



- oo < Vqm < 1. (41) 
The distribution of the quantum potential is given by 

P{Vqm) = ^ 

where 



exp{—iaVQM)0{O')da, (42) 



0ia) = (exp(myQM)) = / d*'^/(^) exp(zaVQM), 



(43) 

/(]?) is given by the same formula as (fT8|) . however, with 

vector Jt = {u,v,Ux,Vx,Uy,Vy) with the same correla- 
tors as (HH). 

For (|43p we may now write with e — 1, which is the 
only case accessible in closed analytic form. 



(a) = — / dudvTxTy exp 
2n J 



, (44) 



with 



^x ^ - duxdvx exp 

TT 



2 _ 2 , iajuUx - VVx) 



(45) 



The same expression holds for Fj^ . The integration in ([l5 
gives 



r r 

x^ y 



-ip 



a — ip 

Substituting (gHI) into (gll) we obtain 

drr^ 



(46) 



eia) 



exp(ia-r72) (47) 



where r — y/p. Finally, substituting that into (|42|) we 
obtain the distribution function for the quantum poten- 
tial 



P{Vqm) = 



1 



2(3/2 - Vqm) 



(48) 



The distribution (gS]) is normahzed as J^^ P{V)dV = 
1. The distribution of P{Vqm) is shown in Fig. [3] and 
compared to a numerical computation of the same statis- 
tics based on the Berry conjecture for chaotic wave func- 
tions in 



ip(r) 



(49) 
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FIG. 3; (color online) The distribution of the quantum po- 
tential (|48[) for e — 1 compared to numerical histogram based 
on the Berry function in Eq. (|49p . 



Here A is the area of the random monochromatic plane 
wave field with |k„p — 1 and the amplitudes for the 
random plane waves obey the relation (a^) = jj. The 
Berry function in Eq. corresponds to e = 1. 



billiard referred to in the introduction, the mass m should 
be the effective conduction band mass m*). Below we 
consider the specific case of small wave lengths A as shown 
in Fig. m We will also comment on the case when A is 
large compared to the dimensions of the cavity. 
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FIG. 4: View of the scattering wave function in the open 
Sinai biUiard for the case (A) shown in Fig.[S]for k'^ = 30.878 
in dimensionless units (see text) and for small aspect ratio 
d/L — 2/67 (ratio between the widths of the leads and the 
billiard). The system is asymmetric because the two opposite 
leads are slightly off the middle symmetry line of the nominal 
billiard. Only the lowest channel is open in the leads. 



V. NUMERICAL SIMULATIONS OF 
SCATTERING STATES IN AN OPEN CHAOTIC 
ELECTRON BILLIARD 

A billiard becomes an open one when it is connected 
to external reservoirs, for example, via attached leads. A 
stationary current through the system may be induced 
by applying suitable voltages to the reservoirs (or by a 
microwave power source as in Section lyl]) . Here we con- 
sider hard-walled Sinai-type billiards with two opposite 
normal leads. A first step towards a numerical simu- 
lations of the quantum stress tensor is to find the cor- 
responding scattering states by solving the Schrodinger 
equation — V^V = ^^"0 for the entire system. The numer- 
ical procedure for this is well known. Thus we use the 
finite difference method for the interior of the billiard in 
combination with the Ando boundary condition [lol . l4l| 
for incoming, reflected and transmitted solutions in the 
straight leads. Once a scattering wave function has been 
computed in this way the fraction residing in the cavity 
itself is extracted for the statistical analysis. To ensure 
statistical independence of the real and imaginary parts 
u and V a global phase is removed as discussed in Ref . Q . 
By this step we also find the value of e. The interior wave 
function is then normalized as defined in Section liiil 

For the numerical work it is convenient to make the 
substitution x x/d and y y/d where d is the width 
of the leads. Here we use dimensionless energy k'^ = 
E/Eo,Eo = n^/{2md'^). (In the case of a semiconductor 



To ensure that the scattering wave function complies 
with a complex RGF we consider a small aspect ratio 
d/L as in Fig. 12 (see also Ref. The actual numeri- 

cal size of the Sinai billiard in Fig. [4] is chosen as: height 
346 (along transport), width (L) 670, radius 87, and 20 
for the number of grid points across the wave guides (d) . 
Within this configuration we now only excite scattering 
wave function with characteristic wave lengths A <C i. 
As expected from Fig.[3]the wave function statistics show 
that both real and imaginary parts, u and v, obey Gaus- 
sian statistics to a high degree of accuracy. Results for 
transmission T and e are are shown in Fig. [5] 

The corresponding distributions for the QST compo- 
nents are given in Figs. [5] and [7] supplemented by the 
distributions for with the x axis directed along trans- 
port. There is indeed an overall good agreement between 
theory and simulations. However, in the statistics for 
in Fig. [S] one notices a tiny difference at small values of 
jx ■ The reason is that there is a net current at this value 
of e, which is not incorporated in our choice of analytic 
isotropic RGFs. The deviation is, however, much too 
small to have an impact on the statistical analysis pre- 
sented here because the net current is such a tiny fraction 
of the entire current pattern within the cavity. The case 
B with e = implies that the scattering wave function in 
the cavity is real (standing wave with transmission T = 
as seen from Fig. [SJ . Therefore there is no current within 
the cavity. 

The agreement with the analytic results for RGFs and 
the present numerical modeling for billiards of finite size 
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FIG. 5: (color online) The transmission probability T (solid 
line) and e (dashed line) as function of the dimensionless en- 
ergy for the Sinai billiard in Fig. |4] Two open circles 
labeled show the case with maximal e = 0.75 (A) and with 
the minimal e = (B). At most only one channel is open in 
the leads. 
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FIG. 7: (color online) Analytic and numerically simulated 
distributions of the components of the QST for the case B in 
Fig. E] (e = 0). The simulated distribution for Tyy is nearly 
identical to PiT^x) and therefore not shown here. Because 
e vanishes there is not any current within the cavity. (The 
choice of lines in the graphs and units are the same as in Fig.|S] 
Because the close agreement between theory and simulations 
differences are hardly noticeable.) 





FIG. 6: (color online) Analytic and numerically simulated 
distributions of the components of the QST and probability 
density current jx along the transport axis for the case A 
shown in Fig. [5] (e = 0.75). As in Figs. [T] and [2] the diago- 
nal components are measured in terms of their mean values 
while Txy and jx are given in terms of \fl^Tiy) and ^/ijV), 
respectively. Solid lines refer to analytic results for RGFs 
(Section IIIII and Ref. 0]) and histograms to the present nu- 
merical modeling. Because of the close agreement between 
the two cases differences are barely resolved. 



is obviously good in the range of energies explored here. 
In order to smooth fluctuations in the distributions of 
the stress tensor we have averaged over the energy win- 
dow shown in Fig. [5] (without scaling e to 1 in contrast 
to Fig. [11] of Section IVip . In this way one finds a per- 
fect agreement between theory and numerical simula- 
tions as shown in Fig. [51 For future reference we note 



that the presence of net currents through the billiard ap- 
pears to have little or no influence on the distributions 
for the present two-lead configuration and choice of en- 
ergy range. We also note that the present results are not 
sensitive to the position of the leads. For example, we 
have also performed simulations for Sinai billiards with 
one dent only and with the leads attached at corners. 





FIG. 8: (color online) Analytic and numerically simulated 
distributions of the components of stress tensor T^x and Txy 
averaged over the energy window given in Fig. [S] The theo- 
retical curves are obtained also by averaging over computed 
e- values shown in Fig. [5] (The choice of lines and units are 
the same as in Fig. |S] The agreement between theory and 
simulations is, however, excellent, hence any small differences 
are not resolved on the scale shown here.) 



We now turn to the complementary case of long wave 
lengths (low energies) . The low energy regime is achieved 
for large aspect ratio d/L which selects wave functions 
with A a few times less than L. Moreover a low energy 
incoming wave often excites bouncing modes. Numerics 
for the case fc^ = 12 and large aspect ratio d/L = 1/7 
show that the scattering wave function may be rather dif- 
ferent from a complex RGF. Hence the corresponding dis- 
tributions for individual states deviate appreciably from 
the theoretical RGF predictions in Section HITl However, 
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by averaging over a wide energy window, as above, one 
closes in on theory. In this way one introduces an en- 
semble that, for practical purposes, mimics the random 
Gaussian case. This aspect may be useful in experimen- 
tal circumstances in which the short wave length limit 
might be hard to achieve. 

VI. EXPERIMENTAL STUDIES 




FIG. 9: (color online) Sketch of the microwave billiard. The 
basic size of the billiard is 16 cm x 21cm. The attached leads 
have a width of 3 cm. The central shaded field (10 cm x 14 cm) 
indicates the region where the data have been collected. The 
measurement grid size was 2.5 mm. The gray regions at the 
end of the two leads indicate absorbers to mimic infinitely long 
channels. The crosses indicate the antennas in the system and 
the winding path illustrates how the third probing antenna is 
moved across the billiard during measurements. 

In quasi-2D resonators there is, as outlined in the in- 
troduction, a one-to-one correspondence between the TM 
modes of the electromagnetic field and the wave func- 
tions of the corresponding quantum billiard [l| . The z- 
component of the electromagnetic field corresponds to 
the quantum- mechanical wave function -0, and the wave 
number fc^ = uP' jc? to the quantum- mechanical eigenen- 
ergy, where lo is the angular frequency of the TM mode 
and c the speed of light. In the present study a rectan- 
gular cavity (16 cm x 21cm) with rounded corners has 
been used, with two attached leads with a width of 3 cm. 
Antennas placed in the leads acted as source and drain 
for the microwaves (see Fig. [9]). Two wedge-shaped ob- 
stacles had been attached to two sides of the billiard to 
avoid any bouncing ball structures in the measurement. 
The same system has been used already for the study of 
a number of transport studies [H, [l3| and for the statis- 
tics of nodal domains and vortex distributions [i^. A 
more detailed description of the experimental setup can 
be found in Ref. |43| . The field distribution inside the 
cavity has been obtained via a probe antenna moved on 
a grid with a step size of 2.5 mm. To avoid boundary 
effects, only data from the shaded region (see Fig.lH]) has 
been considered in the analysis. 

The transmission from the source to the probe antenna 
has been measured on the frequency range from 5.5 to 
10 GHz with a step size of 20 MHz, corresponding to wave 
lengths from 3 to 5 cm. The transmission is proportional 




FIG. 10: The figure shows difi^erent quantities obtained from 
the measurement at the frequency v — 8.5 GHz. In (a) the 
intensity of the wave function is shown and in (b) its phase. 
The plot (c) shows the Poynting vector of the system be- 
ing equivalent to the probability current density in quantum 
mechanics. In (d)-(f) different components of the QST are 
shown, namely xy (d), xx (e) and yy component (f). Dark 
areas indicate higher values. 

to the electric field strength, i. e. to the wave function, at 
the position of the probe antenna. This assumes that the 
leak current into the probe antenna may be neglected. 

To check this we compared the experimentally ob- 
tained distribution of wave function intensities p = 
with the modified Porter-Thomas distribution (see 
e.g. Ref. 0) 

p (iVl') = A^exp (-^^2 |v,|2) /o (/Vm' - 1 |V>|') (50) 
where 

^,= ^(^e+^^ and e'^{v')/{u') (51) 

Here e has not been fitted, but was taken directly from 
the experimentally obtained values for (u^) and (w^), 
where we have ensured that (uv) = by applying a 
proper phase rotation as in and commented on in Sec- 
tion |V|. Whenever x^, the weighted squared difference of 
the experimental data and the modified Porter-Thomas 
distribution, was below Xcutoff=l-l, the pattern has been 
selected for the final analysis of the statistics for the QST 
components. 

Since the wave functions are experimentally known, in- 
cluding their phases, the quantum-mechanical probabil- 
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FIG. 11: (color online) Results for the experimental statisti- 
cal distributions for the components of the QST stress tensor 
obtained by a superposition of all experimental data scaled 
to e = 1 as explained in the text. The solid lines correspond 
to the theoretical predictions in Section Hill for e = 1. 



ity density j = Im^*V7/', and the components of the QST 
can be obtained from the measurement. As mentioned, 
distributions of current densities and related quantities 
have already been discussed previously in a number of 
papers (see e. g. Ref. but the QST has not been 

studied experimentally before. As an example Fig. \W\ 
show intensity (a) and the phase (b) of the measured 
field at one frequency, as well as the probability current 
(c) and different components of the stress tensor (d - f). 

The analysis of the data has been performed in dimen- 
sionless coordinates x = fcr. Since u and v are two inde- 
pendent random wave fields we may rescale the imaginary 
part to obtain e values of one, thus mapping the exper- 
imental result to the situation of a completely open bil- 
liard. This step made it easy to superimpose the results 
from many field patterns of different frequencies which 
originally had different e values. For the analysis all wave 
functions passing the test mentioned above have been 
used. Altogether 83 of 225 possible patterns have been 
taken in the analysis. 

Figure [11] shows the distribution of the QST compo- 
nents obtained in this way. In addition the theoretical 
curves are shown as solid lines. From the figure we see 
that there is a good overall agreement between experi- 
ment and theory, but also that non-statistical deviations 
are unmistakable. 

Deviations between experiment and theory had already 
been found by us in the past in an open microwave bil- 
liard, similar to the one used in the present experiment. 



FIG. 12: Plot of the net current as it is defined in Eq. (l52|l . 
The shaded regions are indicating three different regimes of 
net current strength which had been used in the later analysis. 

in the distribution of current components [l^, [3| . For 
the vertical y component a complete agreement between 
experiment and theory was found, but for the horizon- 
tal X component the experimental distribution showed, 
in contrast to theory, a pronounced skewness. The ori- 
gin of this discrepancy was a net current from the left to 
the right due to source and drain in the attached wave 
guides. In a billiard with broken time-reversal symmetry 
without open channels, a complete agreement between 
experiment and theory had been found, corroborating 
the net current hypothesis. 

For a quantitative discussion of the net current we in- 
troduced the normalized net current for each pattern 




where the average is over all positions in the shaded re- 
gion in Fig.O In Fig. [12] the y component of jnot is plot- 
ted versus its x component for each wave function. One 
notices an average net current pointing from left to right, 
with an angle of about twenty degrees in an upward di- 
rection. For the analysis we discriminated between three 
regimes for the strength of the net current. Additionally 
we performed a coordinate transformation such that for 
each pattern the vector of the net current is aligned along 
the positive a;-axis. This rotation has been done for all 
experimental and numerical results in this section. 

In Fig. [T3] the results for the three different regimes of 
net current strengths are shown. For the distributions 
of the XX and the yy component of the QST, a clear de- 
pendence on the net current strength is found, where the 
deviations from theory increase with an increasing net 
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FIG. 13: (color online) Histograms of the QST distributions 
obtained from experimental data. The thick lines corresponds 
to the smallest net currents (see Fig. I12|) . the thin lines to 
intermediate ones and the dashed lines to ones with the largest 
net current. As in Fig. [Tl]the solid lines correspond to the 
theoretical predictions in Section UlTl for e — 1. 
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FIG. 14: (color online) Histograms of the QST components 
obtained from the simulations according to the wave function 
in Eq. (|53|) . As in previous figure, the thick lines correspond 
to low, thin to intermediate and dashed lines to large net 
currents. 



current. T^y is only slightly affected by the net current, 
if at all. In the limit of a tiny net current, all experimen- 
tal distributions approach the theoretical ones. 

To further test the influence of the net current on the 
distributions of the stress tensor, we performed a numer- 
ical simulation with random plane waves. Each random 
wave field was calculated on an area of 500 mm x 500 mm, 
with a grid size of 2.5 mm. The random wave field con- 
sisted of 500 plane waves with random directions and 
amplitudes. The frequency used for the numerics was 
v = 5 GHz. To introduce the net current we first per- 
formed a random superposition of plane waves according 
to Eq. (|49p . and then added a normalized plane wave 
with the wave vector K' pointing in the same direction 
as the net current observed in the experiment, 

4'{r) = l^a'e'K'-- + a„e*'^"-- j (53) 

The strength of the resulting net current was adjusted by 
a prefactor a' . The best agreement between the experi- 
ment and the numerics was found for a' ~ 0.45. To get 
sufficient statistics we averaged over 200 different wave 
functions. Thus a pattern similar to the one shown in 
Fig. [T^lwas obtained with a cloud of dots extending over 
all three regimes of net current considered with its center 
in the central regime. 

Fig. [13] shows the distributions for the QST compo- 
nents for numerical data derived from Eq. (j53p . The 



same three regimes as for the experimental study have 
been used. The results from this type of simulation are 
in good qualitative agreement with the experimental re- 
sults. In particular the deviations from the theory in 
Section IlIII increase monotonously with the net current, 
just as in the experiment. 

An obvious question is why these net current effects 
are unimportant in the simulations for the Sinai billiard 
presented in Section |Vl One may argue that the num- 
ber of independent plane waves entering at a given fre- 
quency is given by the circumference of the billiard di- 
vided by A/2, where A is the wavelength. Also the width 
of each wave guide is of the order of A/2, i. e. the relative 
net current is approximately given by the total widths 
of all openings divided by the circumference of the bil- 
liard. Following this argumentation the net current in the 
experiments amounted to about 10 percent of the total 
current, whereas in the simulations for the Sinai billiard 
it was smaller by a factor of 10; i.e. too small to be of 
any importance in the simulations. 

We have shown that in the limit of small net currents, 
the distributions of QST components obtained from the 
experiment are well described by means of the random 
plane wave model and the analytic distributions in Sec- 
tion |TTT1 On the other hand net currents are unavoidable 
in open systems. As indicated by the simulations for a 
Sinai billiard in Section |Vl the magnitudes and effect on 
the different stress tensor distributions may be sensitive 
to geometry and energy. Hence it remains an open task 



11 



for theory to incorporate net currents in order to allow 
a more realistic comparison with present experimental 
results. 
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